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Q^ Abstract: The concept of weighted asymmetries is revisited for semi- inclusive deep in- 

(— I elastic scattering. We consider the cross section in Fourier space, conjugate to the outgoing 

'~~' hadron's transverse momentum, where convolutions of transverse momentum dependent 

psj parton distribution functions and fragmentation functions become simple products. Indi- 

^ vidual asymmetric terms in the cross section can be projected out by means of a generalized 

Q^ set of weights involving Bessel functions. Advantages of employing these Bessel weights are 

^j| that they suppress (divergent) contributions from high transverse momentum and that soft 

• factors cancel in (Bessel-) weighted asymmetries. Also, the resulting compact expressions 

f^ immediately connect to previous work on evolution equations for transverse momentum 

^"^ dependent parton distribution and fragmentation functions and to quantities accessible in 

lattice QCD. Bessel-weighted asymmetries are thus model independent observables that 

augment the description and our understanding of correlations of spin and momentum in 

nucleon structure. 
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1 Introduction 

In the factorized picture of semi- inclusive processes, where the transverse momentum of 
the detected hadron Ph± is small compared to the photon virtuality Q^, transverse mo- 
mentum dependent (TMD) parton distribution functions (PDFs) characterize the spin and 
momentum structure of the proton [1-7]. At leading twist there are 8 TMD PDFs. They 
can be studied experimentally by analyzing angular modulations in the differential cross 
section, so called spin and azimuthal asymmetries. These modulations are a function of 
the azimuthal angles of the final state hadron momentum about the virtual photon direc- 
tion, as well as that of the target polarization (see e.g., Ref. [8] for a review). TMD PDFs 
enter the SIDIS cross section in momentum space convoluted with transverse momentum 
dependent fragmentation functions (TMD FFs). However, after a two-dimensional Fourier 
transform of the cross section with respect to the transverse hadron momentum Ph± , these 
convolutions become simple products of functions in Fourier 67^-space. The usefulness of 
Fourier-Bessel transforms in studying the factorization as well as the scale dependence of 
transverse momentum dependent cross section has been known for some time [9-15]. In 
this paper we exhibit the structure of the cross section in 67^-space and demonstrate how 
this representation results in model independent observables which are generalizations of 
the conventional weighted asymmetries [6, 7]. Further we explore the impact that these 
observables have in studying the scale dependence of the SIDIS cross section at small to 
moderate transverse momentum where the TMD framework is designed to give a good 
description of the cross section. In particular we study how the so called soft factor cancels 
from these observables. The soft factor [14-19] is an essential element of the cross section 
that emerges in the proofs of TMD factorization [11, 13-15]. It accounts for the collective 
effect of soft momentum gluons not associated with either the distribution or fragmentation 
part of the process and it is shown to be universal in hard processes [17]. Depending on 
the factorization framework, it appears explicitly in the structure functions and thus in the 
factorized cross section (see Refs. [14, 18]), or it is completely absorbed in the definition 
of TMD PDFs and TMD FFs (see Refs. [15, 19]). At tree level (zeroth order in as) the 
soft factor is unity, which explains its absence in the factorization formalism considered for 
example in Ref. [8]. However, for a correct description of the energy scale dependence of 
the cross sections and asymmetries involving TMD PDFs, it is essential to include the soft 
factor. Yet, it is possible to consider observables where the soft factor is indeed absent or 
cancels out, these are precisely the weighted asymmetries. 

1.1 Overview on vi^eighted asymmetries 

The concept of transverse momentum weighted single spin asymmetries (SSA) was proposed 
some time ago in Refs. [6, 7]. Using the technique of weighting enables one to disentangle 
in a model independent way the cross sections and asymmetries in terms of the transverse 
(momentum) moments of TMD PDFs. A comprehensive list of such weights was derived 
in Ref. [7] for semi-inclusive deep inelastic scattering (SIDIS). A prominent example is the 



weighted Sivers asymmetry, obtained from the differential cross section da according to 

.wisin{ij>h-<t>s) _ 

! d\Phx_\\Phx_\d(t)hd(l3swi{\Phi_\) sm{(j)h - (t)s) {da{(l)h,(l)s) - da{(t)hAs + t^)} 
j d\Phi_\d4>h\Phi.\d4>sWQ{\Phi.\) {da{(j)h,4>s) + da{cl}h,<j)s + 'n-)} 

(1.1) 

where the integrations are performed over the observed transverse hadron momentum 
I Ph± 1 5 the hadron azimuthal angle (f)h and the spin direction (ps of the transversely polar- 
ized target, and the weights are wi = \Ph±\/zM, wq = 1- At tree level and leading twist 
the weighted Sivers asymmetry [7] then becomes, 

^UT - ^ (0)a, . „(0)a. x ' ^^'^f 

Ea Sa /l i^) ^1 (^) 

where /^j: , /} and D^ are transverse momentum moments of TMD PDFs and TMD 
FFs, and Ca is the electric charge for a quark of flavor a. As explained in greater detail 
in Section 5, the moments in Eq. (1.2) are undefined without a subtraction prescription 
for the infinite contributions at large transverse momentum. Here, we propose generalized 
weights, Wn oc Jn{\Ph±\}3T) with Jn denoting Bessel functions of the first kind, and where 
Bt (in units (GeV/c)~^) is a free parameter that represents the Fourier conjugate to |-P/i±|. 
For the Sivers asymmetry, wi = 2Ji{\P h±\BT) / zMBt and wq = Jo{\Ph±\BT)- This gives 
rise to the Bessel-weighted Sivers asymmetry, which reads 






where f^,f , /} and D\ are TMD PDFs and TMD FFs Fourier transformed with 
respect to transverse momentum as defined in the next section. In the asymptotic limit 
Bt — )• 0, we recover the conventional weighted asymmetry Eq. (1.2), and the Fourier 
transformed TMD PDFs and FFs can be identified with the moments in that equation. 
An advantage of the generalized weight relates to the asymptotic behavior of TMD PDFs 
(and TMD FFs). We will see that this provides a regularization of the infinite contributions 
at large transverse momentum as long as we keep B"^ non-zero. Moreover, our analysis will 
show that soft factors appearing beyond tree level cancel out of the weighted asymmetry. 
The rest of the manuscript is organized as follows: In Section 2 we write down the 
general form of the SIDIS cross section in the TMD factorization framework and show that 
the convolutions in momentum space appear as products in Fourier space. For simplicity, 
this discussion is presented at tree level. Modifications needed to go beyond tree level are 
discussed in Section 3. Even though our arguments are quite general, for definiteness we 
use the framework of Ji, Ma, Yuan [14, 18], here referred to as the "JMY" framework. 
TMD PDFs at the level of matrix elements will be considered in Section 4. In Section 5 we 
will consider Bessel-weighted asymmetries, focusing on the Sivers asymmetry as an explicit 
example. Further asymmetries at leading twist are listed in Appendix F. We will also 




Figure 1. Kinematics of the SIDIS process, compare Refs. [8, 22]. 



consider x moments of TMD PDFs and introduce a method to study Fourier transformed 
moments in lattice QCD and compare with experiment. Our conclusions are presented in 
Section 7. 

2 The SIDIS cross section in Fourier space at tree level 



2.1 Elements of the SIDIS cross section 

The lepton-hadron cross section of SIDIS i{l)+N{P, S) 
[4, 8, 20, 21] in the notation of Ref. [8] as 



£{l)+h{Ph)+X can be expressed 



da 

dxg dy dip dzh d(j)h \ Ph± \ d\ Ph± \ 
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(2.1) 



where we assume one photon exchange. L^^^ and W^'^ are the leptonic and hadronic tensors 
respectively, and the vector Ph± is the transverse momentum of the produced hadron in 
a frame where the virtual photon and the target are collinear, e.g. in the target rest frame 
or 7*P center of mass frame. It makes an azimuthal angle (ph with the lepton scattering 
plane defined by the momenta of the incoming and the final leptons / and /' (see Figure 1). 
We define q = l — l' , and q^ = —Q^ is the virtuality of the photon, tp is the azimuthal angle 
of /' around the lepton beam axis relative to S±, in DIS kinematics dip ~ dcps [21]. The 
subscript "_|_" denotes transverse projection in the target rest frame while the subscript "t" 
denotes transverse projection in the light-cone frame. We use definitions for the kinematic 
variables and the ratio of longitudinal and transverse photon fiux e as in Ref. [8] , 
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where M is the mass of the target nucleon. We employ the standard light-cone decompo- 
sition of four-vectors w^ = oo'^'nH^ + uj~n't + Wj.. In the 7*P center of mass frame with the 
proton three-momentum pointing in positive z-direction, the nucleon carries no transverse 
momentum, Pt = 0, and x = p'^ / P^ denotes the momentum fraction carried by the quark 
(parton) of momentum p. Further definitions of kinematic variables and details on the 
leptonic and hadronic tensor are given in Appendix A and Ref. [8]. 



At tree-level of the hard photon-quark scattering process, and to leading order in the 
1/Q expansion, the hadronic tensor can be written in factorized form as [5, 8, 23] 

2MTy^^ = Y,^l fcfpTcfKT6'^^\zpT+KT-Ph±)Tv{Hx,PTh^A{z,KTh''} • (2.3) 

The quark-quark correlator [11, 24] in the above equation is defined as 

^.,{p,P,S) ^ J ^^ e'^-' {P,S\i^^{0)U[a]Mb)\P,S) . (2.4) 

In Eq. (2.3) an integration has been performed over the small component p^ of parton mo- 
mentum to obtain a correlator that depends on light-cone fraction x and parton transverse 
momentum prp, that is, 

^ij{x,PT) = dp- ^ij {p, P, S) 



(2^p e- " 'PT-'^{P,S\i^^{0)U[Ct]Mb)\P,S) 



■ (2.5) 
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The state \P, S) represents a nucleon with four- momentum P and spin polarization vector 
S, and quark fields are located at position "0" and "6" in coordinate space. The gauge link 
U[Cb] ensures gauge invariance of the correlator [23, 25]. It corresponds to a path in b space 
which is determined by the color flow in the hard sub-process [26, 27]. We will discuss the 
details of the definition of the correlator and the role of the gauge link l{[Ci,] in Section 4. 
Analogous expressions define the fragmentation correlator Aij(z,P'j^) (see e.g. [8]). 

2.2 Representation in Fourier space 

In this section, we rewrite the SIDIS cross section and its transverse momentum dependent 
components in coordinate &t space, similar as previously done in Ref. [28]. Here however, 
we take advantage of the rotational invariance of TMD PDFs and FFs. 
First we use the representation of the (5-function 

5('\zp^ + Kr- Ph±) = / ^ e^6H.P.+K,-P.x) ^ (2.6) 

along with the following definitions, 



d'^br 



W^^Phx) = / 7T^ e-*^^-^'^^ W^^^ihr) , (2.7) 



I>,,(x,z&t) = j d^p^e'^^'^-PT $i,(x,pr) 



^ e-P+fe- (P,5|V',(0)W[C,]V^,(6)|P,5) 
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A,, {zM)= I d^Kr e*''^-^^ \j {z,Kt), (2.9) 



to re-write the leading term in the hadronic tensor, Eq. (2.3), in Fourier space 

2MVF^^ = ^ e^ Tr (^{x, zbT)-i^K{z, brh") ■ (2-10) 

a 

The advantage of the bx space representation is clear: the hadronic tensor is no longer 
a convolution of pj. and Kt dependent functions but a simple product of by-dependent 
functions. This motivates us to re- write the entire cross section in terms of the Fourier 
transform 

d^ = f Ih^e-ibr-P'^J ^ y—(l + ^\L Wf^A 

dx^dydi^dzhd^h\Ph±\d\Ph±\ J (2vr)2 \x^yQ^l-e)\ 2xJ ^" j 

(2.11) 
Next, we decompose the correlators <1> and A into TMD PDFs and FFs in Fourier space. 
Using the trace notation (see also Eqs. (A. 8) and (A. 9) in the appendix) 

^Pl = -Tr(^r) , (2.12) 

and restricting ourselves to leading twist projections, we obtain the following structures 
for $ 

$[^^](x,6t) = h{x,bl)-ie'.,"hTpSr,M~ft^^\x,bl), 
|.[^+^'](:e,6t) = SL~giL{x,bl)+ibT-STM~g^^{x,bl), 
^[io'^+i'-] (^^ 5^) = s^ l^ (^^ 52^) + ^ Sl b^M hf}}^ (x, &^) 

+ 1 (b^b!^ + hlg^^'^M^STphfi'\x,bl) 

-ie^%TpMh^^^\x,bl), (2.13) 

where a = 1,2 and p = 1,2. Similarly, we obtain the following structures for A 

A[^"1(z, &t) = Di{z, &t) - i 4^hpShTa zMhDf^'\x, b^) , 
A[^"^'1(z,6t) = ShLGiL{z,bl)-ibT-ShTzMhG^fl{z,bl), 
A[-°~^'1(^,&t) = S^THiiz^b^T) -iShLb''zMhHt!:^\z,bl) 

+ \ {b'^bP, + ]^blg^Py''MlShTpHtT^\^,bl) (2.14) 

~ie^%TpZMhH^^^\z,b'^j.). (2.15) 

For future applications, we have written down the latter decomposition for the more general 
case of a spin-^ hadron; the expression for a spinless hadron is obtained by setting Sh = 0. 
The above decompositions can be deduced from the existing expressions for $ and A in 
momentum space [5, 29], or starting from the symmetry properties of the correlators $ 
and A and a parameterization in terms of Lorentz-invariant amplitudes, see also Section 
4 and Appendix C. The functions /i(x,&y), giL{x,brp), . . . are the Fourier transforms of 



the usual TMD PDFs /i(x,p|,), giL{x,p\), . . .. For a generic TMD PDF called / and a 
generic TMD FF called D, this Fourier transform is given by 

= 2tt j d\pT\\pT\ Jo(|6t||PtI) fix.pl) , (2.16) 

D{z,bl) = f (fKre'^'T-^T D{z,kI) = 2tt j d\KT\\KT\ Jo{\bT\\KT\) D{z,kI) . 

(2.17) 

Additionally, in Eqs. (2.13) and (2.15) not only Fourier transformed TMD PDFs and TMD 
FFs, but also their b^^'derivatives appear, which we denote as 

= ^^ld\PT\\PT\ (^-^y Jn{\bT\\PT\) fi^.Pl) , (2.18) 

b^-){z,bl)^n\{-^^d,.)y b{zMT) 



27r n! 



(z2M^y^'^^"^^'vWrJ '^-^\^t\\^t\) D{z,Kl) . (2.19) 

The functions /, Z), /(") and Z)(") are real valued and /(°) = /, L)(°) = Z). Taking the 
"asymptotic limit" Ibyl — )• on the right hand side of Eqs. (2.19), we formally obtain the 
conventional moments of the TMD PDFs and TMD FFs, /'"•'(x) and D'"'(z) respectively, 

/(")(^,0) = j d'pr (^^J fix,pl) = /W(x) , 

^W(z,0) = 1^2X^(2^) D{x,K'r)^D(-K^)- (2.20) 

Thus we find that the derivatives in by-space are directly related to moments of TMD 
PDFs and FFs. Finally we re- write the SIDIS cross section of Ref. [8] in the 7*P center 
of mass frame with the proton three-momentum pointing in the negative z-direction (so 
called Trento conventions [22]), as 



da 

dxg dy dtps dzh dcph \ Phi. \ d\ Phi. \ 
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^'^UT,L 
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+ V2e(l + e) Sin(20;, - 05) J2(|bT||i'h±|)-^T^'^''"'^"^ 



+ |5,|Ae 
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(2.21) 



The structure of the cross section is what one gets from a multipole expansion in b^- 
space followed by a Fourier transform, see Appendix B. Each of the structure functions 
■^XY z ^"^ ^T-space corresponds to the Hankel (or Fourier-Bessel) transform of the corre- 
sponding structure function F'-^y z ^'^ ^^^ usual momentum space representation of the cross 
section. The combinations si\i{n(j)h + . • ■)Jn{\bT\\Ph±\) and cos(n0/i + . . .) Jn(|^T||-P/i±|) 
act as basis functions of the combined transform to (|P/i_L|, 0/i)-space. Due to the fact 
that the multipole expansion of the physical cross section terminates, only a finite number 
of terms appear in the cross section, with J3 being the Bessel function of highest order. 
The structures J^xYZ ^^^ functions of |&t|, x and z, but no longer depend on the angular 
variables. Introducing a short-hand notation for products 



V[f(-)D(-^)] ^ x^yel{zM\bT\r {zM^lbrir f<^\x,z^bl) D-("'\z,b^T) , (2-22) 



the leading twist tree level analysis in Eqs. (2.10), (2.13) and (2.15) reveals that the Fourier 
transformed structures in the cross section are simple products of TMD PDFs and TMD 
FFs 

F(0) n(0)i 



Tuv.T = P[/r^ D^] , (2.23) 

UT,T 



j^nip-c^s) ^ -r[f^^^) Df»] , (2.24) 

j^^s{<p,-^^) ^ ^j-(i) ^(0)j ^ ^2.26) 

^in(20,) ^ ^^j^Mi) ^±{l)j ^ (2.29) 

_^sin(30ft-</,s) ^ lp[/i^^2) ^±(l)j^ (2.30) 

For completeness, we also list the above results in terms of the momentum-space structure 
functions F-^'y^ of Ref. [8] in Appendix D. Note that TMD evolution equations are typically 
derived in &T-space and are thus obtained in terms of the same (derivatives of) Fourier 
transformed TMD PDFs and TMD FFs that appear in the equations above, see, e.g., Ref. 
[28], where a similar representation of the structure functions in Fourier space has been 
employed. 

3 Beyond tree level 

The formalism becomes more involved once diagrams beyond leading order in Og are taken 
into account. Various strategies have been proposed to address extra divergences that 
appear at the one loop level and higher order [15-19, 30-34]. The development of these 
frameworks for transverse momentum dependent factorization and the establishing of the 
corresponding factorization theorems is an active field of research (see e.g., Refs. [15, 35]). 
The proposed strategies require the introduction of new variables that act as regularization 
scales, and most importantly as it pertains to the content of this paper, the so called soft 
factors coming from soft-gluon radiation. As stated in the introduction, depending on 
the framework, the soft factors appear explicitly in the structure functions [14, 18], or are 
absorbed into the definition of TMD PDFs and TMD FFs (see e.g., Refs. [15, 19]). We wih 
present general arguments that soft factors cancel in weighted asymmetries, independent 
of the specific factorization framework; however for definiteness we work with the JMY 
framework [14, 18], which is based on the ideas of Collins, Soper, and Sterman for the 
factorization of e~^e~ and Drell Yan scattering [13, 30]. Again we consider the structure 
function giving rise to the Sivers asymmetry, 

(3.1) 



The first term in the following referred to as the "TMD expression", dominates in the 
region where \Ph±\ is small, \Ph±\/z ~ Qt ^ Q- The second term is necessary to properly 
describe the structure function for large transverse momentum, where Qt ~ Q, and where 
fixed order perturbation theory and collinear factorization apply. Here H^^l^ is the 

hard part, and S^^' is a soft factor appearing explicitly in the structure function within 
the JMY formalism. It is the same in all the structure functions ^XY z^ ^^^ ^&i. [28]. All 
other structure functions of Eqs. (2.23) - (2.30) need to be modified analogous to Eq. (3.1). 

The term y^™^^'' {Q"^ ,hrp) represents contributions that are relevant only in the 
region of large transverse momentum |-P/i_l| [19, 36]. Qualitatively, this corresponds to the 
very small br region, z\hT\ ^ 1/Q- Since our aim is to study TMD PDFs, we want to 
focus on the region \Phi.\/ z <C Q where we expect them to give the dominant contribution 
if 2;|6r| S> 1/Q. Nevertheless, since we are considering weighted integrals of structure 
functions, the integrals do include the region of very large |-Ph±|. As a result, the Y term 
in Eq. (3.1) is non-zero even if 2:|br| 3> 1/Q. We note that the Y term is expected to be 
particularly important in the case of a "mismatch" between the tail of the TMD term and 
the I -P/i_L I -behavior obtained from the collinear formalism in the regime of intermediate 
|-P/i±|; i-6., M <^ Qt ^ Q- Matches and mismatches between the collinear and TMD 
formalism have been discussed in detail in Ref. [37]. An important example for the case 
of a mismatch is the cos{2(pfi) asymmetry. One possibility to avoid the discussion of the 
y-term is to explicitly cut off the | Ph± \ integrals at some upper value Atmd • This cutoff 
introduces an error in our extracted TMD expression, for which we give an estimate in 
Appendix G.3. Another option is to simply ignore the Y term. This amounts to keeping 
the TMD term in the large |-Ph±| region, i.e., to include the large- |-P/i_l| -tail generated by 
the TMD term, which would otherwise be corrected by the Y term. In Appendix G.3, 
we show that in the z\bT\ ^ 1/Q region of interest this produces an error that falls off 
at least as a fractional inverse power with increasing |&r|- It should be mentioned that 
this estimate of the behavior of the error applies to the Bessel weighting which we discuss 
below. By contrast, no such error estimate exists for conventional weighting with powers 
of l-P/ixl since such integrals are divergent. Better error estimates, or equivalently, a better 
determination of the TMD region in Bt, can be obtained by an explicit treatment of the 
Y term, which we will leave for future analyses. 

In summary, we find that weighted integrals based on the TMD expression alone are 
valid only in a limited range of Bt- Finally, beyond tree level, the product notation 
V[fD] defined in Eq. (2.22) has to be updated to include further dependences on the 
renormalization and cutoff parameters fj,^, p, C ^-^id ( appearing in the JMY formalism 
discussed in more detail below ^: 

r[f(^)D(n^)]^^^Y,el{zM\bT\rizM|,\bT\rr^''\x,z^b^T,^J'^c,p)D''^^^ 

a 

(3.2) 
^The framework of, e.g., Ref. [19], would require analogous modifications within this formalism. 
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4 TMD PDFs at the level of matrix elements 

Apart from introducing the parameters (, C and p the purpose of this section is to review the 
formahsm of Lorentz-invariant amphtudes underlying the decomposition of $ Eq. (2.13). 
In the framework of JMY, the TMD correlator $ itself involves a soft factor S^^' as already 
encountered above, i.e., Eqs. (2.5) and (2.8) need to be modified. In the following, we label 
the unmodified correlators with the subscript "unsub" : 



U[Cb 
(2vr)4 2 



^uLbip,P^s■,v,^,)= I ^e'p-''Up,s\ m m^v]u[^v,b]ri^{b) \P,S) 



^Ssub(^^>5;^,M) 



(4.1) 



The gauge link U[Cb] is essentially given by two parallel straight Wilson lines running out 
to infinity in the direction given by the four-vector v and back again. The definition of a 
straight Wilson line between two points a and b is 



U[a,b] = r exp (-igj^ d^ A^i^] , (4.2) 



where A^{(^) = T^A'^{(^), c = 1..8 is the (matrix valued) gauge field. A transverse link 
connecting these parallel Wilson lines at infinity can be omitted in the covariant gauge 
used by JMY. In case of SIDIS, the direction v = [^",^^,0] is slightly off the light-cone 
direction n_, while for the Drell-Yan process v is slightly off the light cone direction — n_. 
The shift away from the light cone is time-like in the JMY framework and specified in a 
Lorentz-invariant way by the parameter (^, defined by ^^ = (2P • v)'^/v'^. The parameter (^ 
represents a rapidity cutoff parameter [30]. The above correlator can be parameterized in 
terms of real-valued Lorentz-invariant amplitudes. Here we restrict ourselves to the case 
r = 7^^. Reference [29] lists the following structures 

2^unsub - ^ ^2 +P ^3 ^ J^^ ^vPa^ pAi2 + (y.p) 1 



v-F v-F 



The amplitudes B)^ only appear when the dependence of the correlator on the direction v 
is explicitly taken into account, and were not listed in earlier works [1, 5]. Since v represents 
only a direction, the structures above should remain invariant under re-scaling of f , i.e., 
under the substitution v ^ r]v, for any positive real number r/. This has been ensured by 
dividing by powers of v-P in the expression above. The amplitudes can depend on p"^, p-P, 
v-p/{v-P), v'^ /{v-P)'^ = C^^ and the sign of v-P, which we denote with the superscript (+) 
(SIDIS case). For the Drell-Yan process, v-P has the opposite sign (— ). 
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For our discussion below, we make use of a similar decomposition as in Eq. (4.3). How- 
ever, instead of parameterizing the p-dependent correlator $[", we directly parameterize 
the 6-dependent matrix elements ^^^' of Eq. (2.8) in terms of complex-valued amplitudes 
A\ and B^ that depend on b"^, b-P, v-b/{v-P) and (~'^. This parameterization in Fourier- 
space has already been employed in [38, 39] ^. As shown in appendix C, we can deduce 
this parameterization from Eq. (4.3) using the substitution rule p — t- —iM'^b : 

vP v-P 

~{b-S)e>^^^Pp,b^vp 5(+) - .^^{vS)e>'^^Pp,b^vpB\^^ . (4.4) 

In order to connect to the framework of TMD PDFs, we integrate the correlator ^ over 
the (suppressed) momentum component p~ . The integration with respect to p~ reduces 
the Fourier transform with respect to b'^ to the evaluation of <1> at 6"'' = 0. Moreover, in 
the formalism of JMY, the defining correlator of TMD PDFs needs to be modified with a 
soft factor. The modified, ^"-integrated correlator reads 

^™(---^'^'^^C,.)^/;|ye--- /l^e— 

x\{P,S\ mmb]Ti^{b) \P,S) /5W(fo|,^^p) 



€lAb^P^s■v,^/ 



(4.5) 



where xP"*" = p"*" . The soft factor is given as 

5(+)(b^,//^p) = ^ {<d\Tic{U[-oov + b^,bA_\U[b^,b^ + oov] U[oov,G]U[{),-oov] } |0) 

(4.6) 
and involves another time- like direction v = {v^ ,v^ , 0) slightly off the light-cone direction 
n+, controlled by the parameter p = ^Jv^v^ /v^v^ . Note that p^ + 2 + /9~^ = A{v-v)'^ /v'^v^ 
is a Lorentz-invariant expression. Here, the superscript (+) specifies the sign of v-v^ which 
is different for the SIDIS and the Drell-Yan process. 

In the formalism of JMY, the definition of the soft factor S^'^> above both applies to 
the occurrence of S^^> in the TMD PDF correlator Eq. (4.5) and in the structure function 
Eq. (3.1). In the following, we will consider the case F = 7+. The correlator <!)(+) " 1 can 
be decomposed into contributions from two distinct TMD PDFs: 

cl>w^^"](x,p^,p,g,M^C,p) = /l(x,p^;/.^C,p)- '^^^^^^' /l-^T(x,P^;/^^C,P)• (4.7) 



^ In Refs. [38, 39], a difTerent convention for the position of the quark fields in the Fourier transformed 
correlator 4> has been used. These references introduce $ as $(/, P, S,C) = | (-Pi S\ ^(l) U^[Ci] T -0(0) \P, S). 
In Eq. (4.1) we stick to the more common convention of an operator 7/;(0) . . .-(/'(fc). From translation 
invariance follows that the variable b corresponds to —I in Refs. [38, 39]. In particular, our amplitudes 
Ai{b^,b-P, . . .) correspond to Ai{l^, -IP, . . .) of Refs. [38, 39]. 
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Strictly speaking, f^ should also carry the superscript (+) since it has a different sign for 
Drell-Yan measurements [40]. We now use 



bl = -b^ 



b+=o ' P+ 



, o, M2 v+ / 4M2 ^ ^ 



and insert the parameterization Eq. (4.4) into Eq. (4.5). Comparing with Eq. (4.7) allows 
us to write the TMD PDFs /i and /^ as 



/,(.,p|;M^C,p)=2/^e-(^■^)j[°^^Jo(y^^) 



I(+) (h2 up {b-P)R{e) ^-2 ,.2 



X 



5(+)(-62,/.2^p) 

(4.9) 



2 . „2 ^ .^ _ , ,^2 5 /■ d{b-P) ^„(,.p) r ^iz^' 

47r 

(+) (h-i h P (fe-^)fi(C') A-2 ,,2 



/iV(.,p|;.^C.rt = 4M^g^ j ^ ."<-> y^ ^.o(V-i.V 



X 



M+)(-62,A^^p) 



2m2 / ^i^ e"(^-^) r ^(-fe')'^iW-^'P^) 



(2^) ^0 47r ./I^ 



T 



X^^ 



J(+) /"^j^ .p (b-P)K(c2) .-2 .^2 



5(+)(-62,/.2,p) 



(4.10) 



where 



4b = 4^^ + i?(c')i?i+\ (4.11) 

^2]? = A[f - R{e)Bi-'K (4.12) 

We observe that the amplitudes Bi give rise to structures in Eqs. (4.9) and (4.10) that are 
suppressed by their explicit C-dependence as C — >• co, i-e., in the limit of light-like v. The 
structures also disappear in the limit of vanishing nucleon mass M^ — )• 0. Notice that the 
two independent Fourier transforms in each of Eqs. (4.9) and (4.10) naturally connect the 
TMD PDFs to a manifestly Lorentz-invariant framework and reveal x -f-)- b-P and p^ -f-)- b^ 
to be pairs of conjugate variables. 
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The relations of the amphtudes A^^ and A\^2b ^° ^^^ foT-Fourier-transfornied TMD 
PDFs defined in Eqs. (2.17) and (2.19) are given by 

/f (x,6^;/i2,C,p) = 

l,p)J (27r) "" ^2B ^ bT,b^, ^2 'C ,/X^, 



S(+)(6^,/i^ 



(4.13) 



/{!^^^^(x,&f.;/i^C,p) 



-2 ^^^(6:^) ^.,.P) ^(.) f-b^,6.p,(^-^)f5),c-,.- 



5(+)(b^,/i2,p)y (2vr) 



M2 



(4.14) 



We note that the soft factor would need to remain in the integrand if it were also dependent 
on v-b/yv'^, i.e., the "angle" between the Wilson lines and the vector b separating the quark 
fields in the operator. The above result can also be obtained by comparison to the correlator 
$[t^1(x,6t) in Eq. (2.13), where 



l>(+)[^](x,bT,/.^C,p) 



d{b-P) 

(27r)P+ 



J,x{b-P) unsub 



Ki.Ab,p,s■v,^,'') 



5(+)(-62,/i2,/5) 



(4.15) 



5 Bessel-weighted asymmetries 



As stated earlier, transverse momentum weighted asymmetries [5-7] provide a means to 
disentangle the convolutions in the cross section in a model independent way. Generally, 
the conventional weighted asymmetries are given by 



A 



w 

XY 



Jd\P 



h±\ 



'h±\ 



^#sW(|P, 



h±\ 



daxY 



fd\Ph±\ \Ph±\ d(t)hdct)s daxY 
f d\Ph±\\Ph±\d(t)hd^sy^{\Ph±\Ah,(l)s) (daly 



: for XY = UU 



daj^y 



J d\ Ph± I I Ph± I d(j)h d<l)s ( dal^Y + dcr^ 



'XY 



else, 



(5.1) 



where the labels X,Y represent the polarization, "un" (U), longitudinally (L) and trans- 
versely (T) of the beam and target, respectively. The angles (j)s and (ph specify the direc- 
tions of the hadron spin polarization and the transverse hadron momentum respectively, 
relative to the lepton scattering plane. In case of single or double spin asymmetries da^y 
denotes the cross section with one of the polarizations opposite than for da-^y^ such that 
the relevant structure function is projected out from Eq. (2.21). We have introduced the 
short-hand notation W which is a function containing various powers and Ph± as well as 
angular dependences of the form sm{m(j)h =t n(j)s) or cos{m(ph i n<j)s). For the conventional 
weighted Sivers asymmetry, W = wi su\{(j)h — (f^s)-, where wi = \Phi\l zM as in Eq. (1.1). 
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Based on the expansion of the SIDIS cross section in terms of Bessel functions J„ of 
transverse momentum and impact parameter in Eq. (2.21), we exploit the orthogonality to 
generalize the weighting procedure. Now the weighting is of the form 



A^y{St) 



Id\P^ 



h±\ 



h± 



^d</.5W(|P, 



h±\ 



]Bt) daxY 



!d\Ph±\\Phi.\ 
! d\Phi.\\Phi]d^ 



d(j)hd(l)s JQ{\Ph±\ST) daxY 
>hd(l)sW{\Ph±\,4>h,(l)s;BT) (da^Y -dcTxY 



for XY = UU 



Jd\P 



h±\ 



h± 



J<t>s M\Ph±\l3T) [da'xY + dcr^xY 



.t 



: else, 



(5.2) 



where the weight function W corresponds to that of conventional weighted asymmetries, 
except that we replace 



h± 



-^ JJ\P. 



..i^.)n.(Ay 



(5.3) 



As mentioned earlier, taking the asymptotic form of the Bessel function the conventional 
weights [6, 7] which are oc |-P/i±|" appear as the leading term of the Taylor expansion 
of the right hand side of Eq. (5.3). Furthermore we note that the parameter Bt > 
regularizes UV divergences in moments of TMD PDFs and FFs. More importantly, we will 
show that the parameter Bt > allows us to scan TMD PDFs and TMD FFs in Fourier 
space. In fact, the form of Eq. (5.2) already indicates that the weighting implements a 
Fourier-decomposition of the cross section in transverse momentum space. 

Now we summarize the cancellation of the soft factor. We will illustrate this for the 
Sivers Bessel- weighted asymmetry (for details see Appendix F). One can see from Eq. (2.21) 
that the appropriate weight for the Sivers asymmetry is 



W 



2Jii\Ph±\B T) 
zMBt 



sm 



4>s), i.e., wi 



2Ji{\Ph±\BT) 
zMBt 



(5.4) 



corresponding to \Ph±\/zM in the limit |P/i_l| ^ l/0r- Then the Bessel-weighted Sivers 
asymmetry is 






sin(0h-</>s) 



{Bt) = 

jd\Ph^\\Phi_\ 



2Ji(|Pfex|gT) 

zMBt 



sm{(ph 



05) {da^ 



d(j^) 



! d\Phi.\\Phi.\d(t>hd<Ps M\Phi.\BT) {da^ + da^) 
where the axially symmetric denominator is given by 



(5.5) 



2(x' 



XgVQ^ (1-e) 



1 + 



7 
2x, 



d\Ph±\ \Ph±\d<Phdcl)s M\Ph±\BT) 
d\bT\ 



(2vr) 



bT\Jo{\bT\\Ph±\)J'uu,T, (5.6) 
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and from Eq. (2.21) the numerator is 



2a2 ^^/ tM f 2,h{\Pni_\BT) 



2- 



——-\bT\Jl[\bT\\i^h±\)J-UT ■ (5.7) 

Finally, making use of the closure relation of the Bessel function (see Appendix E) we 
obtain for fixed x,y,z, cancellation of the soft factor S~^^^' {^"^ , p) in Eq. (3.1) from the 
Bessel-weighted Sivers asymmetry, 

2Ji(|p^|Br) , . 

^UT,T i^TJ - 

Y.^elH^^^^{Q^p^p)ff^\x,z^Bl-p^,C,p)Df^\z,Bl;p^lp) 

(5.8) 

Some comments are in order. First, if \Bt\ is large enough, the estimate Y^^j!^ ^ ~ 

— 1/2 

Brp can be applied and may indicate that the y-terms are sufficiently suppressed to 
be neglected for practical purposes (see Appendix G.3), which is what we have done in 
the above equation. The above result for the Sivers asymmetry can be generalized to any 
other asymmetry in the SIDIS cross section, Eq. (2.21). We summarize those results with 
the full kinematic dependences in Appendix F. Weighting with Bessel functions at various 
values of Bt thus allows us to map out, ratios of Fourier-transformed TMD PDFs as well 
as azimuthal and spin asymmetries. 

Secondly, the hard scattering factor H^^j.^ ^ is expected to be the same as the 
unpolarized one Hjjjjj,^ because the Sivers effect concerns unpolarized quarks which leads 
to unpolarized scattering on the partonic level. This expectation is confirmed in a recent 
calculation by Kang, Xiao and Yuan [41] at the one loop level, but should hold to all orders. 
Since this feature of the Sivers asymmetry is not shared by the other asymmetries, we will 
stick to writing H^^l^ to avoid potential mistakes. 

Thirdly, it is important to note that in the limit Bt -^ 0, the cancellation of the 
soft factor becomes trivial, since the soft factor S'^{bT, p, p) is unity at bx = 0. This 
has been shown in Ref. [30], but it can also be seen easily from its formal definition in 
terms of Wilson lines given in Eq. (4.6). Using S~^^^\0, p"^, p) = 1 shows that the Phi- 
integrated cross section does not depend on the soft factor, as expected because the collinear 
factorization result should in principle be retrieved (after a proper regularization, which is 
a highly nontrivial matter as discussed in [31]). Due to the asymptotic properties of Bessel 
functions for small arguments, we recover conventional weighted asymmetries in the limit 
Bt^O 



\Ph±\ 



sin((/>h-0s) , 



UT,T y-^Bi'^^yj ~ 

Jlgel Hf<p-^^\Q\p\p) ft^'^\x-p\Cp) dT\z;p^Cp) 
Ea^l H^u,AQ'^P''P) ff\x;p^,Cp) Df\z;p^,(,p) 
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where f-^^ji , /^ , and Z)^ are moments of TMD PDFs and fragmentation functions as 
defined in Eq. (2.20). We caution the reader that these moments are not weh-defined, since 
the corresponding integrals are known to faU off too slowly at large transverse momentum 
[37]. Furthermore, the arguments made earlier that the y-terms are small are no longer 
applicable. 

Lastly, we briefly address what is known about the energy scale dependence of the 
conventional weighted asymmetries. The current knowledge on this is limited to the one- 
loop level. Choosing the factorization scale fi = Q removes the Q dependence from the hard 
scattering function H that is a function of In Q^//i^. This will lead to a Q dependence in the 
transverse moments of the TMD PDFs only [15, 19]. The scale dependence of /} (x; Q^) 
is known, assuming that a proper definition of the TMD PDF can be used, such that the 
zeroth moment corresponds to the collinear function /i(x;Q'^) after the regularization is 
removed. The same applies to D^ (z; Q'^)- For the first moment of the Sivers function one 
can exploit that it is directly related to the Qiu-Sterman function Tp{x, x) [42] as shown in 
Ref. [23]. The evolution equation of the Qiu-Sterman function has recently been obtained 
[43-46] allowing for evolution of the weighted Sivers asymmetry. The evolution of Tp(x, x) 
is not autonomous, since it depends not just on Tf(x,x) itself. This is true even in the 
large- A^c limit, but in the large-x limit it does become autonomous [46, 47]. It indicates 
that /-i^ji (x) evolves logarithmically with Q^ just like /i(x), only falling off faster at a 
given x value as Q^ increases. The evolution has also been calculated for moments of other 
TMD PDFs such as h^ [44, 48, 49] and is similar to that of f^j^ but simpler since 
nonsinglet. In addition, the evolution of the first moment of the Collins function, H^ is 
calculated in [50, 51]. 

6 Average transverse momentum shift and Bessel-weighted counterpart 



In a similar manner to Section 5 we now consider the soft factor cancellation in the average 
transverse momentum shift of unpolarized quarks in a transversely polarized nucleon for a 
given longitudinal momentum fraction x. This shift is considered in [52] and defined by a 
ratio of the pj.- weighted correlator: 



{Py(.x)) 



TU 



f£pTPy^^+^^^^\x,PT,P,S,fl^C,p) 



rfi'T^\x;n'^,C,p) 



5±=0,St=(1,0) a (^!^ '^'^) 



(6.1) 



where f^j^ and f^ are the moments defined in Eqs. (2.20). Obviously, the average 
momentum shift is very similar in structure to the weighted asymmetry Eq. (5.9). While 
the weighted asymmetries are accessible directly from the P/i_L-weighted cross section, the 
average transverse momentum shifts are obtained from the p-p-weighted correlator and 
could in principle be accessible from weighted jet asymmetries. As already mentioned, the 
integrals defining the moments of TMD PDFs on the right hand side of the above equation 
are divergent without suitable regularization. In the following, we therefore generalize the 
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above quantity, weighting with Bessel functions of |py| instead. In particular, we replace 

Py = \Pt\ sin((/)p) — > ^ sin{(t)p - (ps) , (6.2) 

where (/>s = for the choice St = (l^O) in Eq. (6.1). The correlator $(+)[''' ] reads in 
terms of the amplitudes A\ and B^ , 

$(+)[^+l(x,Pr,P,5,/x2,C,p) = jf I'^lbTl { Jo(|foT||PTl)2lS/^ 

- Mlforl \St\ sin(0p - <Ps) Jiilbrl \Pt\) 2^2^/^ } > (6-3) 



where we abbreviate 

(27r) ■ ^ ' 

The Bessel- weighted analog of Eq. (6.1) is thus 

B, _ I d\pT\\PT\ I d<Pp '''^^^f ^h mi<Pp- cl>s) <^^+^^''^Kx,PT,P,S,f,^C,p) 



{Pyi^Wu = 



Jd\pT\ \Pt\ J d<f>pM\pT\BT) <I>(+)[7+](x,Py,P,5,//2,C,p) 

'? ^12B (^~"T'"'-^' M2 )C ;A^ 



|St| = 1 



-M 



fAt^{-Bl,h.P^±I^^X-',,- 



/f)(x,e2;^^C,/^) 

Again, the soft factors cancel. At this point, the independence of the soft factor on v-h/vv"^ 
is crucial. In the limit Bt — )• 0, we recover equation (6.1), {py{x))j^u = {py{x))TU i which 
we have thus shown to be formally free of any soft factor contribution. However, we caution 
the reader again that the expressions at Bt = can be ill-defined without an additional 
regularization step. 

We can go one step further and form ratios that are also integrated in x, with weights 
exp(— zx^Sl). For Bl = 0, this is the same as taking the lowest 2;-moment that appears 
in the Burkardt sum rule [52]. The reason it is interesting to look at such quantities 
is their renormalization properties. Another motivation to discuss such quantities here 
is lattice QCD. Taking x-moments is a standard ingredient in lattice computations of 
nucleon structure, see e.g., Ref. [53] for a review. First exploratory studies of TMD PDFs 
on the lattice [38, 39] focus to a large degree on computations of the lowest x-moment of 
distributions, but access to finite values of Bl is also possible. By "integration over x" we 
mean an integration over the entire support of the correlator; this includes contributions 
from negative x which correspond to anti-quark contributions, see e.g., Ref. [5, 39] for 
details. In particular, the x-integrals of the two TMD PDFs /i and f^ can be decomposed 
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as 



Jo 

(6.6) 
where /i and f^ are anti-quark distributions. In analogy to Eq. (6.5), one can consider 



/dlPrllPrl /#P 2Ji(|p:r|gT) g-^(^^ _ ,^5) J-dxe-*^^^$(+"^+l(x,p^,p,^ 



\Py)TU 



5,M^C,p) 



M 



/cixe-"g^/,y^(x,gj.;/x^Cp) 
/dxe--«i/f)(x,^2.^2^^^p) 

M^:: ^^ - (6 7) 



In this case, the cancehation of the soft factor occurs even for a soft factor that has a 
dependence on v-hjyv^. In the last line of Eq. (6.7), the amplitudes in the numerator 
and denominator parameterize the same matrix element ^^ 1(6, P, 5; f , /x), hence they in- 
volve the same bi-local quark-quark operator. The work on renormalization properties 
of non-local operators involving Wilson lines in Refs. [54-59] suggests that the operator 
'i/'(0)Z//[Cb]'(/'(6) might renormalize multiplicatively for not too small |bT|) compare also [34]. 
As a result, the quantity {py)'j^ ^ would be renormalization scheme and scale independent 
(up to the evolution with the rapidity cutoff parameter C,"^), since all multiplicative renor- 
malization factors would cancel in the ratio. This observation was already made in Ref. [39] 
and is consistent with TMD factorization, which also involves only multiplicative renor- 
malization for l-Pftj.! <C Q or Bt S> 1/Q. For smaller Bt, mixing with gluonic operators is 
expected, as it is known that the Qiu-Sterman function for quarks and gluons (related to 
/it (^) b^ = 0)) mix under changes of the scale [43-46], thereby preventing the cancel- 
lation of multiplicative factors in the ratio considered here. The properties of quantities 
like {Py)j^ ^ remain to be studied more thoroughly. They could be interesting objects to 
make contact between theory predictions from, e.g., lattice QCD, and experiment. 

7 Conclusions 

We have shown that rewriting the SIDIS cross-section in coordinate space displays the 
important feature that structure functions become simple products of Fourier transformed 
TMD PDFs and FFs, or derivatives thereof. The angular structure of the cross section 
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naturally suggests weighting with Bessel functions in order to project out these Fourier- 
Bessel transformed distributions, which serve as well-defined replacements of the transverse 
moments entering conventional weighted asymmetries. In addition, Bessel-weighted asym- 
metries provide a unique opportunity to study nucleon structure in a model independent 
way due to the absence of the soft factor S~^^^' , which as we have shown cancels from these 
observables. This cancellation is based on the fact that the soft factor is flavor blind in 
hard processes, and it depends only on 6-^,//^,/?. Moreover, evolution equations for the 
distributions are typically calculated in terms of the (derivatives of) Fourier transformed 
TMD PDFs and FFs. As a result the study of the scale dependence of Bessel-weighted 
asymmetries should prove more straightforward. For the above stated reasons we propose 
Bessel-weighted asymmetries as clean observables to study the scale dependence of TMD 
PDFs and FFs at existing (HERMES, COMPASS, JLab) and future facilities (Electron Ion 
Collider, JLab 12 GeV). Our results are also easily generalized to other processes where 
TMD factorization is valid, such as e~^e~ annihilation and Drell-Yan processes. 
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A Conventions and useful relations 

In Section 2 we use definitions for the kinematic variables and the ratio of longitudinal and 
transverse photon fiux e as in Ref. [8], 

Q' _ P-q _ P-Ph _2Mx l-j/-^7V , . 

2P-q P-l p-q Q i _ y + ^y2 _^ i ^2^2 

where M is the mass of the target nucleon. The off-collinearity of the process is character- 
ized by the variable Qt introduced through 



qT = q + {l+ ql/Q^)xP - Ph/z, Qt = ^J -q^ (A.2) 

and for Qt <C Q, one finds \Phi.\ ~ zQt-, see, e.g. [37]. The leptonic tensor is 

L^,u = 2{l^C + U' - {l ■ l')g^u + iXee^ua^rq^) , (A.3) 
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where we neglected the lepton mass. The hadronic tensor is 

2MW^^'=Y,fj0^5^'^^ (q+P-Px-Ph) {P, S\ MO) \Px, Ph) {Px, Ph\ MO) \P, S) . 

(A.4) 
For an arbitrary four- vector cj, we introduce the usual light-cone decomposition as 

u}'' =uj^n1+uj-n'^+uj!^ , (A. 5) 

where w = {uj^ it uj^)/^/2 and where the basis vectors n*^ and n^ are 

K = 4^ (1, 0, 0, 1) , n^- = ^ (1, 0, 0, -1) , (A.6) 

such that n± ■ nzp = 1, n± ■ n± = 0, n± ■ lot = 0. Note that cot ■ ^T = —^t- In the 7*P 
center of mass frame with the proton three-momentum pointing in positive z-direction, we 
can decompose the proton and parton momenta as 

p^ = xP+n^ + ^^^n^ + p't , (A.7) 

where x = p'^ / P^ is the quark light-cone momentum fraction. 

Finally, using Eq. (2.12) in Section 2.2, we write the quark-quark correlator recon- 
structed from the trace projections 






1 ^5|,[7''7^] _ 1^^^^^5|,[^<x"^^] + 1;^^[1] ^ (A.^ 



4 "'^ ' 2 " 

where a = 1, 2 and /3 = — , 1, 2. Here we have also included twist-3 and twist-4 terms. For 
the fragmentation correlator we have the following expression, 

- ^7/37'At^'^^] - lza„,7'A[-"'^^l + ilAl^l . (A.9) 

B Multipole expansion and Fourier transform 

This appendix shows the simple underlying mathematical structure of Eq. (2.21). Let us 
treat all kinematic variables except for Phi. and (ph as constants. Consider the cross section 
a in Fourier space as some arbitrary function that depends on by. This dependence can 
be formulated in coordinate space (|br|,</'6), 

oo 

~a{\hT\Ab)= Y. ^''"^'~^n{\hT\) (B.l) 



21 



which is nothing but a multipole expansion with | bx \ dependent coefficients o"„ . Performing 
a Fourier-transform of a back to momentum space (|Ph_L|, (ph), we obtain 

r2TT 



27r /n 27r ■^^ 

•^'^ n=-oo 

^ ,-0, I ^ i^^i (_i)-j„(|P,^||b^|) ^„(|6^|) . (B.2) 



n=—oo 



Using e*"'^'' = cos{n4>h) + isin(n(/)/i) and J^ = (— 1)"J_„, it is evident that the last hne 
of the above equation has exactly the form of the cross section Eq. (2.21), where a finite 
number of the a±n is given by simple linear combinations of the structure functions J'xyz^ 
and the rest is zero. In our case, the Bessel function J„ with the highest n is J3, which 
appears in combination with the angular sin(3(/>/j — (ps) modulation and turns out to be 
associated with the quadrupole deformation of parton densities h^j,. 

C Parameterization of the correlator in 6-space 

First, we briefly review the relevant properties of the correlator under symmetry transfor- 
mations. Applying Lorentz transformations (L), parity transformation (P), time-reversal 
(T) and hermitian conjugation (f) to the matrix elements, we find that the correlator fulfills 

(L) ■■ ^^2sn^{p,P,S;v,f^) = <l>i'li^^^/^'(Ap,AP,AS;A^,/.) , (C.l) 

(P) ■■ ^Ssub(P> P^ S; V, ^i) = ^^22X^ {p, P, -S; V, /i) , (C.2) 



.[r] 



hW'r*r'-/^. 



(T) ■■ [KUiP, P, S; V, ^l)\ = d^V^s^b "'iP, P, S; -v, n) , (C.3) 

(t) ■■ [€lAp^P^S■v,^,)\* = <i>l7;r/l(p,P,5;^;,/x) . (C.4) 

where we denote the sign change of spatial components of a given vector c; that is, c = 
{(P , — c^, —(? , — c^). From hermiticity (f) follows that the Ai and B^ in Eq. (4.3) are real 
valued. Time reversal (T) does not constrain the number of allowed structures, because 
it changes the sign of vP. Instead, time reversal (T) establishes relations between SIDIS 
amplitudes A\ , B^ and Drell- Yan amplitudes A\ , B^ ' . 

For any of the transformations T G {L, P, T, f}, the Eqs. (C.1)-(C.4) are of the general 
form 

nmp,w)) = '^{Tp{p),TUw)) (C.5) 

where we have omitted the subscript "unsub" and the renormalization scale /U, and where 
the symbol w summarizes all dependences on F, P, S and v. Here 7$ is either the identity 
function or complex conjugation. The transformation rule 7p{p) maps onto Ap, p or p 
and thus fulfills a-b = Tp{a)-Tp{b) for any two vectors a and b. The Fourier-transformed 
correlator 

$(6, w)= I (fp e-'P-^^{p, w) (C.6) 
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transforms according to 



$('^r,(6),r.H). (C.7) 



For example, $ transforms under hermitian conjugation as 



;[r] 



[^ort^o] 



(t) : [KL^{h,P,S-v)\ = $;7_J K-b,P,S-,v) . (C.8) 

Let f{p,w) be any of the structures preceding the invariant amphtudes in the parame- 
terization of $. The structure f{p, w) is a homogeneous function of some degree n in 
p, i.e., f{ap,w) = a^f{p,w) for any number a. For example, the structure f{p,w) = 
Ai(v-P) (P'S)^'^'^"'^ ^i^Pa'^/s preceding Bg in Eq. (4.3) has degree n = 2. If we define 
f (b,w) = f{-iM^b,w), then 

n [f{b,w)) =n{-iM^rnifib,w)) = f {n{-iM^)Tpib),%,iw)) = f n^ . 

(C.9) 

This shows that / transforms like $ in Eq. (C.7). We conclude that the parameterization 
of $ can be found by the substitution p — )• —iM'^b in the structures parameterizing <I>, and 
we arrive at Eq. (4.4). The amplitudes A^ and B\ introduced this way are no longer 
constrained to be real valued functions. Instead, hermitian conjugation Eq. (C.8) yields 
the relation 



r(+)/'/,2 A. p ,,.?,//'„. p\ /'-2 ,,2^ 



(+)r^2 _h.T3 _o,.hl(o,.-p\ ^-2 ,,2^ 



A^^^h^ b-P, vb/{v.P),C', m')J = A^>{b', -b-P, -vb/{vP), r^ /x^) . (C.IO) 

D Structure functions in terms of Fourier transformed TMD PDFs and 
FFs 

The structure functions of Ref. [8] can be expressed in terms of Fourier-transformed TMD 
PDFs and FFs as 

Fuu,T = x,Y.el Jj^\bT\ M\bT\ \Ph±\) /r(x, z'bl) DUz, b^) , (D.l) 

p^<p-^s) ^ -a:,J2^l I ^\br\' M\br\ |P.x|) Mz ~ftT^'\x,z'bl) bl{z,bl) , 

(D.2) 
FLL = x,Y.el j '^\bT\M\bT\\Phi.\)riL{x,z''bl) DiizMr) , (D-3) 
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(D.4) 

fJ^(^'^+^-) = x,5]e^|^|6Tp Ji(|foT||P.J) M,z hUx,z'b'r) Ht^^'\z,b'r) , 

(D.5) 
FuT"''^ = x,^ell^\bTfJ2{\bT\ \Phi.\)MM,z' ht^^'\x,z'b'^) Ht''''\z,bl) 

(D.6) 
F^'''^^ = x,Y.-lj^^\^T?U\bT\\PHi.\)MM,z^h^^^^^^^ 

(D.7) 

j^s\n{'i<f,h-<t>s) Sr^ 2 [ 'A^t\,, \i t nu I ID l\^^^'»^^I-La{2)/ 2. 2 n £V-La(l)/ . 2 \ 

^C/T =XB2^e^J T^l^^l H\bT\\Ph±\) ^ \t (a^'^ ^t) ^1 (2,^t) 



,2^ 



(Di 



E Cancellation of the soft factor in the Sivers asymmetry 

Making use of the closure relation of the Bessel function 

r°° 1 

J d\Ph±\\Ph±\Jn{\Ph±\\bT\)Jn{\Ph±\BT) = -^d{\bT\-BT), (E.l) 

we obtain for the expression in Eq. (5.6) 

I d\Ph±\\Ph±\d(l)hd(t)sM\Ph±\BT) I ^\bT\M\bT\\Ph±\)J'uu,T 
= x^^el Huu,TiQ^f^^p) jd\Ph±\ \Ph±\ jd(Ph jd<PsJoi\Ph±\BT) 

X / ^ \bT\ M\Ph±\ IbrDfT^ix, z%h /^^ C, P) S^-'Hbh ^^ p) bf^\z, bh P, C, P) 
= 27rx^ Y,el Huu,TiQ^,P^p) f^^^x.z^B^; t,\C, p)~S^^\Bh^^\ P)Df^"{z,Bl; p,Lp) 

a 

(E.2) 
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Next, we consider the following expression in the numerator of the asymmetry, Eq. (5.7), 
I d\PH±\\P,±\ I d<p, I dcPs ^^'^l^'^^^^^ sin^(0. - cPs) 

xMzft^'^''{x,z'bht,^Cp) S^^\hl,p\p) Df^\z,hl,pMp) 
= 27rx,J2elHi^^%--''^\Q\p^p)ft^'^^{x,z'Blp^C,p) 

a 

X SM ^Bl p\ p)Df'^\z, Bl p\ C/z, p), 



(E.3) 



where we have used the closure relation Eq.(E.l), and 



2tt i-2tt 

cos^ {m (ph + n (j)s)d(f)s = sin^ {m (ph + n(j)s)d(ps = t^ , (E.4) 

'0 Jo 

for integer n and m. Thus, we obtain 

^^ (l-e) 1^^"^ 2xgJ 

EaelH^u,TiQ'^P'^P)fT''i^^^^'3^T-'P'^C^P)S^iBhp'^p)D?"iz,B^^-p^t^ ' 

(E.5) 

F Bessel-weighted asymmetries 

Here we introduce the Bessel weights Wn 

2 ^" 



Wn = Jn{\Ph±\l3T)n\{ — ] , (F.l) 



and summarize the Bessel-weighted asymmetries at leading twist: 
Double Spin 



Jo(|Phx|BT)/K ^ _ o?^ (1-^) V 2xg 



^^^'i + ^K/r^^ 



yQ2 (i_e) l^i -r 2x^ 

Y.aelH,,{Q\p\p)~gfl\x,z^Bl-p\Cp)D^?\z,Bl-p\Lp) 
j:^elH^^^{Q^,p^,p)ff^''{x,z^B^^;p^X,p)Df''''{z,B^^;p^,C,p) 



(F.2) 
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Worm Gear 



^Q^ (l-e) l^^ "^ 2a;g 

(F.3) 
Collins 






S/Q2 (i-e) 



1 + 



7 



^^ Eg e^ ggy^^^"^ (Q^ ^^ P) /^f ^"(x, z^gf.; /i^ C, p) nt'^'^^z, Bl- p\ C, p) 
Ea el Huu,TiQ', I^^P) fT"i^, z^Bl- p\ C, P) Df^\z, B^; p?, (, p) 

(F.4) 

Boer — Mulders 






j,Q2 (i-£) ' 1 + 

^^ Eg el HZ^^^-\q\ p^p) ht^'^^jx, z^Bl- p\ C, P) Ht^'^^jz, B^; p\ (, p) 

Ea el Huu,TiQ', /"^ P) /f ^"(x, z^B^^; p^ (, p) of^^z, B?,; p^ (, p) 

(F.5) 

Kotzinian — Mulders 






e 



AuL '' " (^t) = 2 



_a£_ y^ (11 7- 
yQ2 (1-e) \^^ 2xg 

., Eg el Hfjf'^-\Q\ p\p) /.jj^)"(x, Z2g2 ■ ^2^ ^^ ^) ff^±(l)a(,^ ^2 . ^2^ ^ ^) 

Ea e^ ^c/i/,t(Q^ m'> p) fT"(^, ^'^t; ^^ C, P) Df^\z, Bl; p\ (, p) 

(F.6) 

Pretzelosity 



^C/T V^rj - ^ 1 



J/Q2 (i„£) ^^i -^ 2a;^ 



^^ Eg e^ ggf ^'--'^^^ (Q^ M^ p) /.j-f "(X, Z2g2 ■ ^2^ ^^ p) ^^±(l)a(,^ ^2 ■ ^2^ ^^ ^) 



F.7) 
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G Suppression of high transverse momentum 

G.l Suppression of the tail in Bessel-weighted integrals 

In this section, we are interested in the contribution from the high-momentum region to 
Bessel-weighted integrals. We begin by deriving convergence criteria and upper bounds for 
an integral of the form 

poo f'OO 

In,Ai^,A^)= duOjJn{ojOA{oj)=C^ du U Jni^) A{u / ^ ■ (G.l) 

In the subsections to follow, uj will assume the role of a momentum, Ipj^l, \Kx\ or |P/i_l|, 
while E, will represent |for| or Bt- In the equation above, A is a placeholder for a given 
function of (jJ. We restrict our discussion of the integral /n,A(C) ^w) to the region A^^^ S> 1, 
where the Bessel function in the integrand can be approximated by 



•^nM~y — sin(z^ + J-^) for iy»l 
First, let us consider a function A that fulfills the condition 

|A(a;)| < ca;^° for any w > A^^ 



(G.2) 



(G.3) 



where c > and a are two real valued constants. Using the envelope of the Bessel function 
l'^n(j^)| ^ vW^i '^6 obtain 



4,A(e,Aa 



-2 






OO 

OO 

A„? 



dv 



vJn{y)A{v/i) 



dv\l — c — 



vr v- 



a 



'2A| 



cA^" ~ e 



-1/2 



(G.4) 



where the last equal sign only holds if the integral is convergent, i.e., if a > 3/2. 

Next, consider a function B for which ^/ujB{uj) is monotonously falling in the region 
UJ > A(^ and converging to zero for a; — )■ oo. In this case, we can make use of the oscillatory 
behavior of the Bessel function to show convergence of the integral In,B{i^^w)- Let N 
denote the smallest possible integer such that A^^^ + 7r/4 — n7r/2 < Ntt. We decompose the 
integration according to 



/■N-K 

In,BitK) = r^ / 

^A„5 



NTT-TT/4+mT/2 



dvuJn{y)B{v/i) 



Tl 

°^ r{N+j+l)n-n/4+mT/2 

+ r'y^/ duu.ui^)B{i,/o 

j^Q J{N+j)TT-7T/4+mr/2 

^ V ' 

T2 



(G.5) 
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Applying the mean value theorem, we can find v G [A^j^, A^vr — 7r/4 + 27r/2] such that 

r-N-K—-K/i+n-K/2 

K C'^V9B{v/i)\j- ( cos(A^g + 7r/4 - ra^/4)^ - cos(A7r) ) (G.6) 

from which we derive bounds for the first term: 

< (-l)^+iri < 2^C^^B{v/0 < 2y|r VAc.^ + 2^^(Ac.) . (G.7) 

Using the mean value theorem again to determine the points Vj, the second term becomes 

T2 ^ C"^ Y] y^jB{i^j/0 / du J- sin(zy + 7r/4 - mr/2) 

•^0 J{N+j)iT-TT/4.+nn/2 V VT 



00 



= 2V ^r' E ^/i^B{^?,/0{-l)'''-' ■ (G.8) 

The expression on the right is an alternating series that fulfills the Leibnitz-test for con- 
vergence and is bounded by the size of its first term, 

<(-l)'^r2 < 2^^r^V^B{uo/0< 2y^rVAc.e + 2^S(A^). (G.9) 

Since Ti and T2 have opposite signs, we arrive at a combined upper bound 



/n,B(e,Aa.)| < 2y^rVA..e + 2^5(A^) « 2W^5(a^) ~ r^/\ (g.io) 

In summary, we find that the integral Eq. (G.l) converges for any function A{u}) that 
decays faster than w"^' ^ which in turn determines an upper bound of the integral of order 
^~^'^, Eq. (G.4). The requirement for convergence can be relaxed to functions decaying 
faster than w^^" if monotony of ^/uj times the function is ensured. In this case, Eq. (G.IO) 
gives an estimate of an upper bound that decays with ^^^". We remind the reader that 
these bounds are only valid for ^ ^ A^^^, but for all n. 

G.2 Fourier-transformed TMD PDFs, TMD FFs and their derivatives 

Using the mathematical results from the previous sub-section, we investigate which of the 
(derivatives of) Fourier-transformed TMD PDFs /("^(x, b^) and TMD FFs D^"-\x, b^) are 
well-defined by the right hand sides of Eq. (2.19). Their behavior in the high transverse 
momentum region has been studied in detail in Ref. [37]. They find power-suppressed tails 
of the form 

f{x, \px\) ~ -j 1 — X "logarithmic modifications" , (G.ll) 
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for integer powers m. Analogous expressions hold for the TMD FFs. Comparing the right 
hand side of Eq. (2.19) with the criterion for functions of type B, we find that convergence 
is maintained if n < m — 1/2. The logarithmic modifications do not play a significant role 
since logarithms grow more slowly than any polynomial. 

The analysis of Ref. [37] reveals that (up to logarithmic modifications) fi,giL, ^i, /"*", 
9i,hT,h^,fT,9T,hL,h,eL,e,fL,9-^,eT,e^,Di,D-^,G-^,H,E ~ l/p'^. For these func- 
tions, the corresponding zero-derivative and single-derivative Fourier-transforms f^'^\x, brp), 
f^^'{x, brp), D^^\z, brp) and D^'^'{z, bp) exist. A second group of distributions exhibits the 
high- momentum behavior f^, gix, hj^p, h^ , fj^ , gp , h^p^, H^ ~ ^/Pt- Fo'^ these latter func- 
tions, the existence of n-derivative Fourier-transforms /'"•'(x, b^) and D^'^'[z, by) is ensured 
up to n = 3. Again, we point out that these results are only valid for {bxl > 0, while the 
limiting case I&tI =0 leads to divergent integrals [37]. 

G.3 Systematic errors from the region at large Ph± 

TMD frameworks have been designed to give a good description of the cross section at 
low transverse momentum, i.e., for \Ph±\/z ^ Q. However, in weighted asymmetries we 
integrate over the whole range of |-P/i±|. The contributions from high |P/i±| thus lead to 
theoretical errors in the results if one does not have a description of the cross section that 
is valid there, even when one restricts to the region z\bT\ ^ 1/Q- The Y term can in 
principle be included to eliminate those errors, but its Fourier transform is expected to be 
power suppressed in the region z\bT\ ^ 1/Q, because it was shown to be power suppressed 
at small |-P/i_l| [13, 36]. Dropping the Y term means that we approximate the full result by 
the large |-P/i_L|-tail of the TMD expression. This in general may be a bad approximation, 
but the question is whether it will affect the result much for z\bT\ ^ 1/Q- In addition, 
extending the integrals to arbitrarily large transverse momenta ignores the fact that the 
physical cross section should vanish above a certain maximum transverse momentum value 
|-P/i±|max (see also Refs. [12, 36]). In this appendix we are going to estimate the effect of 
these various simplifications. 

The Y term will be significant only in a finite region of | Ph± \ '■ between a scale Atmd 
and \Ph±\ma.x- Note that both these scales will depend on Q. We can bound the error from 
neglecting the Y term in terms of its maximal value. As long as |6t| ^ ^tmd ^ l-f /iJ-lmaxi 
we can approximate the Bessel-function as in Eq. (G.2) to obtain, 



Y^-/j<^^'^+-\Q2^b'T)^ J d\P,^\ |p,^| 27rJ^(|foT||P.±|) Y^^/r^^'^-^-^Q^Pl^) 

/"''"'""^''dip,^! |p,^| 2nM\bT\\PH±\) yl'^g°^(^'^'^+■•■)(g^pL) 

. (G.12) 



>-TMD 

' 27r 



^ (|-Ph±|max - Atmd) 2^/ 

V I Or Atmd 



T^sin/cos (TV </>,,+...) 
^XY,Z 



Here \Y^yz [max is the maximum absolute value of Y in the range between Atmd 

and l-P/ixImax- It can be estimated from the (perturbatively calculable) y-term. Thus, 
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I 
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( "TMD I region". Qr<S Q )i^ "Y region", Qt ~ Q 
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Figure 2. Schematic illustration of important scales for Bessel- weighted asymmetries before and 
after the Fourier-transform. 



Eq. (G.12) shows that the theoretical error from neglecting the Y term is (at least) sup- 
pressed as Ibrl"^/^. An explicit treatment of the y-term in Eq. (3.1) could eliminate this 
theoretical error to a given order in a^ in the Fourier transformed TMD PDFs and TMD 
FFs extracted using Bessel weighting. We will not do this here. 

The second error coming from extending the TMD expression beyond | PhL \ max is more 
suppressed and therefore less of a concern. Following a similar procedure as before we can 



-iSin/cos(Af(^h+...)l 



estimate it to be suppressed as |br|~ / . Let [-F^y^"°*'""''''^"''']TMD denote the structure 
functions as determined purely within the TMD framework, i.e., from convolutions of TMD 
PDFs, TMD FFs and a potential soft factor. The contribution to its Fourier transform 
coming from the large | PhL \ region can be bounded using that the TMD expression (times 



'hi. 



1/2N 



is a monotonically decreasing function of |-P?j_l|. Thus, applying Eq. (G.IO), 



|Pfe±|i, 



d\PHi\ \Pki\ 2vrJ^(|6r||P,^|) [Ff^/™^(^^'^+■■•)]TMD(Q^PL) 



< 4, 



•^'^1-' /i±|max 



|6t| 



^sin/cos(Af0h + ...)^ 
XY,Z 



- j-iaiii/i;uH(^iy (/>h-r...Jl /^2 \ T3 |2 \ 



(G.13) 



where the upper bound applies as long as Iforl ^ \P 



hi.\ 



This second error is therefore 



far less important than neglecting the Y term. The reason this same behavior could not 
be obtained for the Y term is that it is not expected to be a monotonically falling function 

of |Ph±|. 

Finally, let us consider what error would be introduced if all | -P/i_l | integrations of the 
experimental data were to be cut off at Atmd- In this case, we would be able to use Eq. 
(G.13) as an error estimate, except that |-P?t±|mix would need to be replaced by Atmd- 
Again the error estimate would be valid provided Ibj-I ^ ^tmd ^'^'^ provided the structure 
function times |-P/i±|^ is monotonically falling, i.e., in its tail region, beyond Atmd- This 
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simple cutoff method is expected to be useful when Q^ is very large, such that Atmd can 
be chosen large with confidence. 

In Fig. 2 we illustrate how the contributions from the TMD and the Y dominated 
regions contribute to the Fourier transform. The contributions from the region |-P/i±| > 
Atmd are only suppressed in the region of large Bt > 1/Atmd- Therefore, an analysis 
without Y term at too low values of Bt has to be considered with caution. However, also 
the region of large Bt has to be treated with care in case of Bessel weighting, as one starts 
to probe the oscillations of the Bessel function. This is relevant whenever 1/Bt becomes 
smaller than the experimental resolution in transverse momentum. A finite transverse 
momentum resolution |Ph±|res can, for example, be a result of binning of the experimental 
data, as indicated in the figure. 
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